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Abstract 



Q I We consider winding number transitions in the two dimensional 0(3) non-linear 

sigma model, modified by a suitable conformal symmetry breaking term. We discuss 
the general properties of the relevant instanton solutions which dominate the transi- 
tion amplitudes at finite energy, and find the solutions numerically. The Euclidean 
period of the solution increases with energy, contrary to the behavior found in the 
abelian Higgs model or simple one dimensional systems. This indicates that there is 
a sharp crossover from instanton dominated tunneling to sphaleron dominated ther- 
mal activation at a certain critical temperature in this model. We argue that the 
electroweak theory in four dimensions should exhibit a similar behavior. 
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Gauge theories of the strong and electroweak interactions are characterized by 
a multiple vacuum structure. Tunneling transitions between different vacua are re- 
sponsible for physically interesting effects, such as baryon number violation in the 
electroweak theory. At zero temperature and energy these winding number transi- 
tions are dominated by the familiar zero energy instanton solutions of the Euclidean 
field equations, with vacuum boundary conditions. 

At finite temperatures thermal activation over the potential barrier separating 
the multiple vacua can occur in addition to quantum tunneling. The static classical 
solution whose energy is equal to the top of this barrier between neighboring vacua is 
the sphaleron solution. At sufficiently high temperatures transitions between different 
winding number sectors are dominated by classical thermal activation with a rate 
controlled by the energy of the sphaleron. 

In simple one dimensional systems there is a smooth crossover from the zero energy 
instanton dominated tunneling transition to the high temperature sphaleron domi- 
nated regime. The corresponding classical solutions interpolating between these two 
situations are known as periodic instantons with turning points at finite Euclidean 
time 13, which give a non-pertubative contribution to the partition function and tran- 
sition rate at temperature f3~^. One would expect similar considerations to apply in 
quantum field theory, although the situation is much less well explored and very few 
classical solutions of this kind are known. 

Periodic instantons appear naturally also in the context of zero temperature field 
theory, if one considers the transition probability between different winding number 
sectors at fixed energy, E. This probability may be expressed in the form, 

P{E)=Y.\{f\SVE\^\\ (1) 

where S is the 5-matrix, Ve is a projector onto energy E and the initial and final 
states \i) and (/| lie in different winding sectors. The path integral for P{E) is 
saturated by a periodic instanton solution [jl|, and the probability is given, with 
exponential accuracy, by the exact analog of the quantum-mechanical formula, 

P(E)~exp[E/3-5(/3)], (2) 

where S{f3) is the Euclidean action of the periodic instanton. The Euclidean period 
P is related to the energy by 

^-- 

The initial and final multi-particle states can be read off from the analytic continua- 
tion of the periodic instanton into Minkowski time at its turning points. Hence peri- 
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odic instanton solutions to the classical Euclidean equations contain non-trivial infor- 
mation about multi-particle transition amplitudes between different winding number 
sectors at finite energy. It has been suggested that by suitably modifying the bound- 
ary conditions in the complex time plane, information about transition amplitudes 
from few particle initial states to many particle final states may be obtained as well 
0. For these reasons it is interesting to study more thoroughly the possible classical 
periodic instanton solutions in field theories with winding number transitions. 

Because of the absence of analytical methods for finding periodic instanton so- 
lutions to the Euclidean equations in interesting field theories, one must rely on a 
numerical approach. In earlier work periodic instantons have been studied numeri- 
cally in the two dimensional abelian Higgs model ^. In this letter we present the 
results of a numerical study of periodic instanton solutions in the two dimensional 
non-linear 0(3) sigma model, which differs from the abelian Higgs model by the ab- 
sence of a scale parameter (before the addition of any 0(3) symmetry breaking). As 
a result, zero energy instantons in the unbroken 0(3) sigma model can have arbitrary 
size, exactly as is the case for unbroken SU{2) gauge theory in four dimensions. As 
we shall see, this leads to a dependence of the period of the periodic instanton on 
energy which is quite different from that observed in systems with only one degree of 
freedom or the abelian Higgs model, and a sharp rather than smooth transition from 
quantum tunneling to thermal activation at a certain critical temperature. 

The two-dimensional 0(3) sigma model has the Euclidean action 

So = ^,Jd'x{d,n,)\ (4) 

where na{x), a = 1, 2, 3 are three components of a unit vector, UaUa = 1. The model 
is well known to possess instanton solutions 0], which may be pictured geometrically 
by first conformally mapping two dimensional Euclidean space R^ onto S"^ with all 
points at infinity identified, and then taking the field configuration which is just the 
identity mapping of this S'^ onto the S'^ of n^ field space. Hence the zero energy 
instanton covers the entire sphere pictured in Fig. 1 as Euclidean space and time 
vary over the infinite range [— oo,+oo]. Because of the conformal invariance of the 
action (^, the zero energy instanton can have arbitrary scale p. This implies that the 
energy barrier between winding number sectors, proportional to 1/p, can be made 
arbitrarily small and no finite energy sphaleron solution exists in the symmetric 0(3) 
model. This is precisely as in the pure Yang-Mills theory in four dimensions. 

In the electroweak theory conformal invariance is broken in the Higgs sector which 
then makes possible the existence of a classical sphaleron solution. In the 0(3) sigma 
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X = p/4 




Figure 1: Geometric representation of the periodic instanton. 



model the conformal invariance can be broken by adding to the action @ the term 



S, = ^Jd'x{l + m) 



(5) 



which also violates the 0(3) symmetry and fixes the vacuum state to be n"'^'""'^^ = 
(0, 0, —1). With the conformal symmetry breaking term (^, the modified model with 
total action S = So + Si possesses the classical sphaleron solution 



ni'^'Kx) 



-2-4 — ^,0,-1 + -^ 
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which has the energy 



E 



sph 
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Geometrically, this solution maps the infinite spatial line onto a great circle beginning 
and ending at the south pole in Fig. 1. 

The spectrum of perturbations around the sphaleron has exactly one normalizable 
negative mode with the eigenfunction Ua{x) = (0, ch~^(ma;), 0) and the eigenvalue 
tu^ = — 3m^. Because of the existence of this negative mode it is possible to determine 
the form of the periodic instanton for energies just below the sphaleron energy (|^). 
Adding a small amplitude e of the negative mode to the static sphaleron gives a 
Euclidean time-dependent periodic instanton solution, 



nt'^^Hx) + e sm(\u-\T)ua(x) 



(8) 



with an energy just below the sphaleron energy. The period of the periodic instanton 
approaches the critical value (3crit = 2n/\LU^\ = 27r/(-\/3m) ^ 3.628/m as e ^ and 
the energy approaches the sphaleron energy, E — > Egph. 

In the other limit of very low energies, E <^ Egph, the periodic instanton solution 
can be found by considering an infinite chain of instanton-anti-instanton pairs on 
the Euclidean time axis. The infinite chain acts as a set of images for the field 
configuration in one link with fundamental period (3. The classical action of this 
approximate solution over the fundamental period is 

The corrections to the action are suppressed by factors of m'^f3'^, m^p^ and p^//?^ 
which all go to zero for E <^ E^ph, as we shall see in a moment. Substituting this 
action into the expression for the rate of tunneling transitions at fixed energy, (§) 
gives 

P(E) ~exp(-^Vr) (10) 



^(/5)^4(l-7r2|^-mVln(mp)) . (9) 
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with 

ly = 87r - g'^(3E - 8n^^ - STcm^p^ In(mp) (11) 



where the values of the period (3 and the instanton size p are determined by the 
saddle-point equations |jl| 

'''' '"" 0. (12) 



dj3 dp 



In the limit E <^ £'„„/, we find 



sph 



27r2 1 1/2 

/3(^) ^ £:^Tr7z:7-77?T ^ (13) 



which implies 



m V\Yi{m/g'^E) 
^^^^ ^ 25/4^ {\n{m/g^E)} ' ^^^^ 

wiE)^8.-'-^\^'^;r, (15) 

^ ^ m \\n{m/g^E)} ' ^ ^ 

and justifies the neglect of the corrections to the action in the limit of very small 
energies. 

This approximate behavior of the periodic instanton at low energy leads us to 
observe that the conformal symmetry breaking term in the action (^) serves only to fix 
the scale of energy and instanton size. It does not distort the periodic instanton from 
the configuration obtained by a simple linear superposition of zero energy instanton- 
anti-instanton pairs of the model with unbroken symmetry. This is possible only 
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because P{E) goes to zero as E ^ 0, i.e., P{E) is an increasing function of energy, 

d(3{E) 



dE 



> . (16) 



Notice that this behavior is nevertheless consistent with the hnear superposition or 
dilute gas approximation because the parameter that controls the validity of this 
approximation is p^//5^, and p goes to zero even faster than P as E ^ 0. Hence it is 
already clear that the behavior of the periodic instantons in the 0(3) model is very 
different from that of systems with only one quantum mechanical degree of freedom 
or the abelian Higgs model. In those models f3{E) —>■ oo as E ^ 0, and the periodic 
instanton contributes to the transition rate as the temperature j3~^ goes smoothly 
to zero 1^. Periodic instantons in the present model with period going to zero as 
E ^ cannot contribute to the low temperature transition rate, although they can 
contribute, and in fact dominate transitions between non-thermal low energy states. 
We shall see presently that the behavior of the period with energy (|l^) persists up to 
the sphaleron energy. 

At energy comparable, but not very close to the sphaleron energy, the periodic 
instanton solution has to be found numerically. For that one must solve the Euclidean 
field equations 

- {dl + dDria + m^Saz - An, = 0, (17) 

together with the constraint, 

riaTia = 1 . (18) 

Here, A is a Lagrange multiplier enforcing the constraint, which is easily found by 
multiplying ([T7|) by n, and using (|^). The periodic instanton solution we seek has 
vanishing time derivative at initial time r = 0, evolves to another turning point at 
half-period (3/2 where it reflects and then returns to its initial configuration at t = (3 
by simply reversing the sign of all r derivatives. Thus we enforce the half-period 
boundary conditions, 

dna{T = Q,x) _dna{T = l3/2,x) _ 

This removes the time translation invariance of the solution. The other boundary 
condition we require is that at spatial infinity the solution approaches the vacuum. 
Since we shall work in a finite box of length L we require. 



na[T, X = -L/2) = n,{T, x = L/2) = n'^^'^^ = (0, 0, -1) (20) 

for rriL ^ 1. 



-m 




Figure 2: The finite rectangular region (r, x) of coordinate space that maps onto the shaded 
region of the sphere in Fig. 1. 

Geometrically, the periodic instanton solution obeying these boundary conditions 
maps the rectangular region of (r, x) coordinate space pictured in Fig. 2 into the 
shaded region of the sphere in Fig. 1. From these figures it should be clear that 
the solution may be chosen to have well-defined symmetry properties under reflection 
through the lines bisecting the rectangular region of Fig. 2, so that at x = 0, nx 
vanishes while ?t,2 and ris reach extrema, and at r = /5/4, n2 vanishes while n\ and ns 
reach extrema. Imposing these symmetries fixes completely the spatial translational 
invariance and rotation invariance around the ns axis of the solution, as well as reduces 
the region we must consider to only one quarter of the full rectangle in Fig. 2. Hence 
we must solve the differential eqs. (|l^) subject to the boundary conditions. 



X = -L/2 
t = 
x = 

t = (3/A 



rii = 0; n2 = 0; 123 = -1; 
drii/dt = dn2/dt = dn^/dt = 0; 
111 = 0; dn2/dx = dn^/dx = 0; 
^2 = 0; drill dt = dn^/dt = 0. 



(21) 



The lattice version of eqs. (|T^) can be obtained by a variational principle starting 
from the discretized action. Let (tj, Xj), where i = 0, . . . ,1 and j = 0, . . . , J , be the 
coordinates of lattice sites. The lattice version of the action (El-H) reads 






(22) 
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where we have used the notations dti = tj+i — tj, dti = {1 /2){dti_i + dti) and similarly 



for dxj, dxj. 



The equations which follow from the action ( p2D can be written in the form 

[S"(r^ - n''n%j = 0, (23) 

where 

B°: = i+^'-^' _ 7-1 -J _ ^''^J+1 _ ^*».J-i _^ m^^as (24) 

Eqs. ( P5D actually contain only two independent equations since its projection onto 
the vector nf^ is identically zero. Together with the constraint equation (p!8D these 
comprise a complete set of three independent equations associated with each lattice 



site. From eqs. (23) it is clear that i?" must be parallel to n" at each point {i,j)- 
Since n"' is normalized to unity, the eqs. (p3D-(p4D can be rewritten in the equivalent, 
symmetric form, 

u^Vb^ + B" = 0, (25) 

for every {i,j). This form is more convenient for numerical calculations. 

To obtain a numerical solution of ( ^5]) we use Newton's method. That is, we 

1. choose an initial field configuration as a first guess; 

2. linearize the equations in the background of the initial field; 

3. solve the linearized equations to obtain an improved configuration; and 

4. iterate until the procedure converges. 

The choice of initial configuration is guided by the known behavior of the periodic 
instanton solution near E = and E = E^ph. For example, for E just below Egph 
an initial configuration of the form (|^) with e ~ 0.5 works well. The convergence of 
the scheme is quadratic, and double precision accuracy is typically reached in 5 to 7 
iterations. 

The solution of the linearized equations in step (3) of the algorithm requires an 
inversion of the matrix of small fluctuations about the trial configuration. Negative 
eigenvalues of this matrix are treated on the same footing as positive eigenvalues, 
and pose no special problem for Newton's method, which is important for the present 
application. Because of the boundary conditions which fix the translational and 
rotational symmetries there are no zero eigenvalues of the matrix near the desired 
solution, which otherwise would be disastrous for this method. 

For a J X J spacetime grid the matrix to be inverted has Pj^ elements. In general, 
it takes of order Pj^ operations to invert such a matrix. However, the sparseness of 
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the matrix makes it more efficient to follow a procedure of forward elimination and 
back substitution along either the x or r directions instead. That is, starting with 
one edge of the region in Fig. 2, such as r = 0, we solve for each r slice of the grid 
in terms of the successive two r slices until the edge r = /3/4 is reached, where the 
boundary condition determines the unknown quantities. Then we reverse direction 
and solve for the unknowns on the previous r slices successively. This allows for the 
matrix to be inverted in order IJ^ operations (or PJ operations if the x direction is 
chosen) , and speeds up the algorithm considerably. In practice, a grid size of order 
100 X 100 can be handled on a typical workstation, which already provides reasonable 
accuracy (~ 1%) for the problem at hand. 

The numerical method described works quite well and yields a periodic instanton 
solution for energies in the range of 1 to 8 in units of mj g^ . Smaller energies require 
finer lattices to obtain the same accuracy, since the instanton size becomes smaller. 
The results of our calculations are summarized in Figs. 3. They were obtained on a 
40 X 40 grid in a spatial box of size L = 8 in units of vrr^ . 

Fig. 2 illustrates a typical periodic instanton solution. As one can see, at the 
turning points r = and r = /9/2 the field configuration is non-vacuum, and the 
components of the vector n" cover the shaded patch of the sphere in Fig. 1. 

Fig. 3a shows the dependence of the action of the periodic instanton, in units of 
1/(7^, on its half-period, as the latter varies from /3/2 ~ 1.3 to /5/2 ~ /5crit/2. The 
action grows smoothly from a value close to 5*0 = 47r to PcritEsph/'^ ~ 14.5. 

Fig. 3b shows the dependence of the half-period 13/2 as a function of energy. Near 
the sphaleron energy (8 in units mj g^\ the period reaches the critical period ficrit- 

Fig. 3c shows the exponent, W = S — Ef3, in the suppression factor for the 
tunneling probability P{E), as a function of energy. W varies from a value slightly 
below 47r at low energies to zero at -E ~ Egp^. 

The behavior of all quantities in these graphs is smooth, so that greater accuracy 
may be obtained in a straightforward way by increasing the number of lattice sites. 

The periodic instanton solutions we have found numerically dominate the winding 
number transition rate at fixed energy. At fixed temperature, the transition rate is 
given to exponential accuracy by 

r(T)= / dEe-^'^PiE), (26) 

JO 

where T is the temperature. Evaluating this integral by the method of steepest 
descent (which is valid in the limit of arbitrarily weak coupling (^^ -^ provided T is 
coupling independent) gives the saddle point condition 

T'l = -^ = (3{E) . (27) 




Figure 3: (a) Action of periodic instanton versus half-period. 



p/2 




Figure 3: (b) Half-period versus energy for the periodic instanton. 




Figure 3: (c) The suppression factor W (in P{E)) as a function of energy. 



However, evaluating the second derivative of the exponent at this saddle point yields 
which tells us that this extremum is a minimum of the exponent rather than a max- 



imum, for the periodic instanton solutions we have found obeying (|T6|). Hence this 
saddle point does not give the dominant contribution to the rate F, and the periodic 
solutions we have found are not the relevant ones for finite temperature transitions. 
If these periodic solutions are the only ones obeying the saddle point periodicity 
condition (|27|) , then the integral in eq. (|2^) is saturated either at i^ = 0, which 
implies 

r(r) ~ exp(-2S'o) = exp(-87r) , 

or at -E = Esph^ where the exponent W{E) reaches zero and the exponential suppres- 
sion in P{E) disappears. In the latter case one finds 

r(T) ~ e^-Esph/T) . 

The crossover between these two cases occurs when the two contributions become 

equal, i.e., at 

rj, _ Esph _ m 
"'' 2So~ n ■ 
At temperatures above Tcr the transition rate is dominated by the thermal activation 
over the top of the barrier, and the sphaleron solution plays a dominant role. At 
temperatures less than T^r, the transition rate coincides with the zero energy tunneling 
rate. At temperatures T < Tcr there would then be no exact classical solution which 
dominates the finite temperature transition rate, and we would expect instead that 
the field configurations which saturate the rate at T < Tcr are the same as those 
saturating the zero temperature transition rate, viz., the constrained instantons 
which are not exact solutions to the Euclidean field equations. Indeed, at T = 
the absence of any exact real finite action solutions immediately follows from scaling 
arguments (Derrick's theorem), and at finite T imposing periodic boundary conditions 
at Euclidean time T~^ does not alter this argument as long as T -C rra. 

The other possibility is that the real periodic instanton solutions we have con- 
sidered in this letter are not the only solutions to the field equations satisfying the 
periodicity condition (pT]), but that there are other solutions, perhaps complex, for 
which the periodicity decreases with increasing energy, and which do dominate the 
finite temperature rate r(r). Since we know that at sufficiently low temperatures 
there are no other real periodic solutions than the ones we have found, any other 
solutions would have to be complex, in order to evade the scaling argument. Indeed 
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it is interesting to speculate that the reason the constrained instanton configurations 
can dominate the zero temperature transition rate even though they are not exact 
solutions is that they are a good approximation (in their action) to some exact com- 
plex solution to the field equations, yet to be discovered. If this speculation is correct, 
then we would have to reconsider the question of whether the real periodic solutions 
obeying the boundary conditions (^) are the correct configurations maximizing the 
transition probability at fixed energy as well. 

We would like to conclude by stressing the analogy between the two dimensional 
sigma model with the symmetry breaking term (^) and the spontaneously broken 
SU{2) gauge theory. In both cases, in the unbroken phase, instanton solutions exist 
and can have arbitrary size. In the broken phase, exact real instanton solutions do 
not exist, and the dominant contribution to the zero temperature (energy) transition 
rate comes from approximate solutions, the constrained instantons. Like in the sigma 
model, in the spontaneously broken SU{2) gauge theory the period of the periodic 
instanton increases with energy at low energies, where perturbation theory is appli- 
cable |l|]. Hence the situation is in all important respects exactly the same as in our 
broken 0(3) sigma model, and one may expect that the finite temperature transition 
rate in the SU{2) + Higgs system is saturated either by zero energy tunneling or by 
sphaleron activation, with a sharp crossover between these two regimes at tempera- 
tures of order Mw, or else new complex periodic solutions to the equations of motion 
of the SU{2)- Higgs system exist which are still awaiting discovery. 
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